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A PROBLEM OF BERRY AND KNOTTED ZEROS 
IN THE EIGENFUNCTIONS OF THE HARMONIC OSCILLATOR 

ALBERTO ENCISO, DAVID HARTLEY, AND DANIEL PERALTA-SALAS 


Abstract. We prove that, given any finite link L in R®, there is a high- 
energy complex-valued eigenfunction of the harmonic oscillator such that its 
nodal set contains a union of connected components diffeomorphic to L. This 
solves a problem of Berry on the existence of knotted zeros in bound states of 
a quantum system. 


1. Introduction 

In [2], Berry conjectures that there should be complex-valued eigenfunctions 
of the harmonic oscillator in whose nodal set ^“^(0) has knotted connected 
components, and raises the question of whether there can be eigenfunctions of a 
quantum system whose nodal set has components with higher order linking, as in 
the case of the Borromean rings, see Fig. Furthermore, Berry remarks that it 
should be possible to construct these sets so that they are structurally stable in the 
sense that any small enough perturbation of the corresponding eigenfunction (in 
the norm with fc ^ 1) still has connected components in the nodal set that are 
diffeomorphic to the knot or link under consideration. 

As a side remark, let us recall uttn] that a physical motivation to study the 
nodal set of a quantum system is that it is the locus of destructive interference of 
the wave function. It is related to the existence of singularities (often called dislo¬ 
cations) of the phase lm(log'0) and of vortices in the current field Im(i/; Vi/>). The 
existence of knotted structures of this type, especially in optics and in fluid mechan¬ 
ics, has recently attracted considerable attention, both from the theoretical HIS] 
and experimental mm viewpoints. 

The main result of this paper solves these problems of Berry by showing that 
any finite link can be realized as a collection of connected components of the nodal 
set of a high-energy eigenfunction of the harmonic oscillator, and that the link is 
structurally stable in the same sense as above. Specifically, we have the following 

Theorem 1.1. Let L be any finite link in K^. Then one can deform it with a 
diffeomorphism $ of so that $(L) is the union of connected components of 
the nodal set where ip is a complex-valued eigenfunction of the harmonic 

oscillator inMf^. Furthermore, the link is structurally stable for the function ip. 

We recall that the eigenfunctions of the harmonic oscillator are the square- 
integrable functions ip satisfying the equation 

— Alp -b \x\'^ip = Xip 


( 1 . 1 ) 
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Figure 1. The problem involves showing that there are high- 
energy eigenfunctions ■i/' of the harmonic oscillator realizing links, 
e.g. the trefoil knot and the Borromean rings depicted above, in 
their nodal set 


in It is well-known that the eigenvalues are of the form 

X = 2N+ 5, 

with N a nonnegative integer, and that the degeneracy of the corresponding eigen- 
space is 5(-/V + 1)(-/V -I- 2). 

The key idea of the proof of Theorem o is that, using techniques introduced 
in UlSli, one can prove that there are complex-valued solutions to the Helmholtz 
equation 


Aip + (fi = 0 


in K^, such that the link L is a union of connected components of the nodal set 
(^“^(0), up to a diffeomorphism. This is pertinent to the study of the eigenval¬ 
ues of the harmonic oscillator because, in balls of radius the high-energy 

asymptotics of the eigenfunctions are determined by the Helmholtz equation (see 
e.g. [IIH]). Heuristically, one can understand why this is true by introducing the 
rescaled variable x := x, in terms of which Eq. (1.1) is read as 


Axtp + 1 /’ = 




The way to make this precise is by computing the high-order asymptotics of the 
Laguerre polynomials, which govern the radial part of the eigenfunctions of the 
harmonic oscillator. Going over the fine details we will see that the accidental 
degeneracy of the eigenvalues of the harmonic oscillator is an essential ingredient of 
the proof too, essentially because it ensures the existence of families of isoenergetic 
eigenfunctions with a rich behavior in the angular variables. 

The proof of Theorem o is given in Section although the proofs of two 
technical lemmas are relegated to Sections and ^ To conclude this paper, in 
Section we will state and discuss a higher-dimensional counterpart of the main 
theorem that can be proved using the same argument. 
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2. Proof of Theorem 11.11 

Let us begin by fixing an orthogonal basis of eigenfunctions associated with the 
harmonic oscillator Hamiltonian. Specifically, we will take 

( 2 . 1 ) ■= , 

where (r, 9, cj)) are spherical coordinates and we are using the standard notation 
for the Laguerre polynomials and the spherical harmonics. Here the indices of the 
eigenfunctions range over the set 

fc^O, 1^0, 

and the eigenvalue corresponding to V’fcim is 

Xki ■— 4/l -p 2/ T 3 . 

Notice that the eigenvalue is independent of m. 

In the following lemma we will describe the behavior of the eigenfunction ipkim 
and its gradient for large values of /c. To state this result, we will use the notation 
Cr '■= xjr for the unit vector in the radial direction and denote by Vg2Y{9,(j}) the 
gradient (in the unit sphere) of a function Y(9, (j)) of the angular variables. The 
proof of the lemma is given in Section 

Lemma 2.1. Let us fix some integers I and m as above. Uniformly for r ^ R, the 
eigenfunction tpkim admits the asymptotic expansion 

ipkim{x) = Aki [ji{\/hdr) + 0(1)] Yim{9,(l )), 

'^f’klm{x) = y/^lAkl [jz'(\/A^T) + 0(4)] Yim{9,(p) Cr 

+ Aki |_Jz( V Xki r) + 0( j:)J - - -, 

as fc —>■ oo. Here ji is the spherical Bessel function of order I and Aki is a nonzero 
constant. 

In the following lemma we construct an even complex-valued solution of the 
Helmholtz equation: 

Aip - 1 - (^ = 0 , 

in such that the link L is a union of connected components of its nodal set up to 
a diffeomorphism. The function ip is smooth and is conveniently given by a finite 
sum of spherical Bessel functions and spherical harmonics. We observe that ip is not 
square-integrable but it has an optimal decay rate at infinity among all solutions 
to the Helmholtz equation. The proof of this lemma is presented in Section]^ and 
exploits ideas introduced in I1I3IS]: 

Lemma 2.2. There are finitely many complex numbers cim, with 0 ^ I ^ Iq and 
—l^m^l, such that the complex-valued function 

Iq i 

^ ^ cim.3i{r)Yim.{9,(j)) 
l—O m——l 

has the following properties: 

(i) The function ip is even, so cim = 0 for all odd 1. 



4 


ALBERTO ENCISO, DAVID HARTLEY, AND DANIEL PERALTA-SALAS 


(ii) There is a diffeomorphism ofM.^ such that $i(L) is a union of connected 
components of the zero set 

(iii) is structurally stable. More precisely, let S be a compact set con¬ 

taining <i)i(L). Then there is some d > 0 such that for any function ip' 
with \\ip — < 5 one can find a diffeomorphism $2 o/K^ such that 

$2 o is a collection of connected components of (^'“^(0) that are 

contained in S. 


Let us take a large integer k that will be fixed later, and which we assume to be 
larger than ^. For each even integer I smaller than 2k we set 

(2.2) ki := k — - , 
so that the eigenvalue 

(2.3) A:=A^^;=4fc + 3 

does not depend on the choice of 1. The desired eigenfunction i/i of the harmonic 
oscillator can then be derived from the function (p constructed in Lemma |2.2| by 
setting 

1=0 m =-1 fed 


for a large enough number k. Notice that, by construction, ip is a, smooth complex¬ 
valued function in L^(IR^) that satisfies the Eq. (1.1) with A as in (2.3). Here we 
have used that cim = 0 for odd I, since the number ki defined in (2.2) is an integer 
only for even 1. 


Let us fix some R > 0 such that the ball centered at the origin and of radius R, 
which we will denote by B, contains the link <I>i(L). We claim that for any <5 > 0 
one can choose k large enough so that 


(2.4) 


iP 




<(5. 


CHB) 


This is a rather straightforward consequence of Lemma |2.1[ Indeed, using this 
lemma, an elementary computation shows that 




V x'f’klr. 


= A 


kl 


= A 


kl 


ji{r)Yim{0,(p) + 0[^ 


j'lir) <^) e, + ^ cf) + 0(i) 


Hence, substituting these asymptotic expressions in the sum for ip we find 

In I 




l—O m— — l 
Iq I 


-^ip.. 


kil 


'\Vx 


- ji(r)Yim{0,4>) 


“ X! X! czmO(i ) ^ 


C 


l—O rn=—l 


c 

k — k 
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provided k is much larger than ^ and \x\ < R. An analogous argument shows 


Va;V'( - ^x^{x] 


C 

k 


so the estimate ( |2.4[ ) follows provided k is large enough. 

ensures that, if <5 is small enough, the function 


Item (iii) in Lemma 


2.2 


has a collection of connected components in its nodal set {ii)[-/'/\) = 0 } given 
by the link ^2 ° where $2 is a diffeomorphism of and $2 o <i)i(L) is 

contained in B. Item (iii) also ensures that the link $2 o $i(L) is structurally 
stable for the eigenfunction ip. This implies that the rescaled link $3 0 <I> 2 od)i(L) is 
a union of structurally stable connected components of ip~^{0), where $3 denotes 
the diffeomorphism of given by the rescaling 

3 ( 2 :) 

The theorem then follows by setting $ $3 o <I >2 o ‘hi- 

Remark 2.3. It is worth noting that the fact that the function ip is even was key 
to constructing the radial quantum number fc; via Eq. (2.2 1 . A straightforward 


modification of the argument enables us to consider the case where ip is odd. As 
is well known, all eigenfunctions of the harmonic oscillator must have a definite 
parity. In particular the nodal set of the eigenfunction ip contains (at least) two 
copies of the link d)(L) as the link d)i(L) constructed in the proof of Lemmais 
contained in the positive octant of K^, which implies that so is ‘h(L). Moreover, 
$(L) is contained in a small ball of radius RX~^/^. 


3. Proof of Lemma 1271] 

The lemma essentially follows from Hilb’s asymptotic formula for the Laguerre 
polynomials [m Theorem 8.22.4]: 


" rT'^’^"(r^) = Afe/jKi/A^r) + 0(fc‘2") 


d 

dr 


J 7-^+5 ! 


e Tr‘Lr^(r 2 )] = [AuJiiV^ir) + 0{k^)\ , 


with 


Am := 


(^) 


-I 


r(fc + ; +1) 


/n \ 'Z J k\ 

This formula holds uniformly for r ^ R. (In fact, the formula for the derivative 
does not appear in the above reference, but it is standard —and easy to prove— 
that this asymptotic formula can be derived term by term). 

The asymptotic expansion for ipMm written in the lemma follows from the iden¬ 


tity ( 2 . 1 ) and the fact that the constant Am can be estimated for large k as 

2 7+1 . 7-1, 

Am = ~~/=^ ^ + Oipk 2 ). 

This is an elementary computation using Stirling’s formula for the factorial and the 
identity 

Y{k + l + ^= + + 

+ 22fe+2J+2 (it + / + 1)! ■ 
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4. Proof of Lemma [ 2 T 2 ] 

Let be a ball centered at the origin that contains the link L. There is no loss 
of generality in assuming that L is contained in the positive octant of that is, 

L C -B n {xi > 0,X2 > 0, X 3 > 0} . 

An easy application of Whitney’s approximation theorem ensures that, by perturb¬ 
ing the link a little if necessary, we can assume that it is a real analytic submanifold 
of 

Let us denote by the connected components of L, with the index a taking 
values in a finite set A. Each component La is an analytic closed curve without self¬ 
intersections. Our next goal is to write the curve La as the transverse intersection 
of two surfaces and 

Since any closed curve in has trivial normal bundle m, there exists an 
analytic submersion 0 q, : Wa —> , where Wa is a tubular neighborhood of La and 

= ^a- We can then take the analytic surfaces E^ := 0“^((—1,1) x {0}) C 
Wa and E^ := 0,^^({O} x (—1,1)) C Wa- Since 0^ is a submersion, these surfaces 
intersect transversally at La = E^ n E^. 

Now that we have expressed the component La as the intersection of two real 
analytic surfaces E;!^ and E^, we can consider the following Cauchy problems, with 

Aui + ui=0, uij^j^=0, 

Here denotes a normal derivative at the corresponding surface. The Cauchy- 
Kowalewski theorem then grants the existence of solutions u^, to this Cauchy prob¬ 
lem in the closure of small neighborhoods 11^ of each surface E^. We can safely 
assume that the tubular neighborhoods Lfa C Ua are small enough so that the neigh¬ 
borhoods corresponding to distinct components are disjoint. Now we take the union 
of these pairwise disjoint tubular neighborhoods, 

y := U (C/I n c/2), 

aeA 

and define a complex-valued function ip on the set U as 

■■= ui + iul. 

The following properties of p are clear from the construction: 

(i) (p satisfies the equation 

Ap + p = 0 

in the tubular neighborhood U of the link L. We can assume without loss 
of generality that U is contained in the positive octant B C {xi > Q,X 2 > 
0, X 3 > 0}, as is the link L. 

(ii) U can be taken small enough so that the nodal set of p is precisely L, i.e., 
L = p-\0). 

(iii) The intersection of the zero sets of the real and imaginary parts of <,3 on L 
is transverse, i.e.. 


(4.1) 


for all X G L. 


rank(V Re(,3(a;), V Im(^(a;)) = 2 
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Let us agree to say that a subset of is symmetric if it is invariant under the 
inversion x i-A —x, and denote hy —U the image of the set U under this map. Since 
U is contained in the positive octant, U n —U = 0. 

Let us then define an even function ip' in the symmetric set 

U' :=UU -U, 


as 


p'(x) := 


p{x) if X G U , 
p{—x) if a; G —U . 

By construction, p' satishes the Helmholtz equation 
(4.2) Ap' + p' = 0 

in U' and its nodal set consists of L and its mirror image under the inversion 


—X. 


Denote by S' a symmetric closed subset of U' whose interior contains the link 
L. Our next goal is to construct a solution of the Helmholtz equation in that 
approximates the local solution p' in the set S. To this end, let us take a smooth 
even function y : —>■ K equal to 1 in a neighborhood of S and identically zero 
outside [/', and define a smooth extension pQ of the function p' to by setting 
Pq •= x^'■: which is an even function too. Denote by 


G{x) 


cos \x\ 
47r|x| 


the Green’s function of the operator A + 1 in which satisfies the distributional 
equation 

AG + G = -(5o 

with So the Dirac measure supported at 0. Since po is compactly supported, we 
obviously have 


(4.3) 


po(x) = I G{x — x)p{x)dx 


with p := — Ai^o ~ ‘^o- The complex-valued function p is even and its support is 
contained in the set U'\S. Therefore, an easy continuity argument ensures that 
one can approximate the integral (4.3) uniformly in the compact set S' by a finite 
Riemann sum of the form 


(4.4) Pi{x) := pjG{x-Xj). 

3 = -J 

Specifically, for any 5 > 0 there is a large integer J, complex numbers pj and points 
Xj G U'\S such that the finite sum ( |4.4[ ) satisfies 

(4.5) llv^i -‘/5'||co(s) < <5. 

By the symmetry of the integrand, these quantities can be chosen such that po = 0, 
P-j = Pj and x-j = —Xj for j > 0, thus guaranteeing that pi is an even function. 
Here we have used that po = p' in S. 

In the following lemma we show how to “sweep” the singularities of the function 
Pi in order to approximate it in the set S by another function p 2 whose singularities 
are contained in the complement of the ball B. The proof is based on a duality 
argument and the Hahn-Banach theorem. 
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Lemma 4.1. For any ^ > 0, there is a finite set of points {zj}j^_j, in 
complex numbers Cj such that the finite linear combination 

(4.6) >P2{x)-- ^ CjG{x-Zj) 

j=-J' 

approximates the function tpi uniformly in S: 


(4.7) 

Moreover, 


\\F2 - i^lllco(S) < ^ ■ 


Co = 0, 


Z-j = —Z-i 


C—j — Ct , 


for all j > 0, so that ip 2 is an even function. 


and 


Proof. Consider the space lA of all complex-valued functions that are linear combi¬ 


nations of the form (4.6 1 , not necessarily even, where Zj can be any point in . 
and the constants Cj take arbitrary complex values. Restricting these functions to 
the set S, U can be regarded as a subspace of the Banach space (7° (S') of continuous 
complex-valued functions on S. 


By the Riesz-Markov theorem, the dual of C°(S) is the space M{S) of the finite 
complex-valued Borel measures on whose support is contained in the set S. Let 
us take any measure p. G M.{S) such that /jja fdp. — 0 for all f G U. Let us now 
define a complex-valued function F G as 


F{x) := ( G{x — x) dfi{x), 

Jr3 

so that F satisfies the equation 


AF + F = -p,. 

Notice that F is identically zero on M.^\B by the definition of the measure p and 
that F satisfies the elliptic equation 


AF + F = Q 


in M^\S, so F is analytic in this set. Hence, since K^\iS' is connected and contains 
the set by analyticity the function F must vanish on the complement of S. 

It then follows that the measure p also annihilates any complex-valued function of 
the form pj G{x — Xj) because, as the points Xj do not belong to S, 

0 = PjF{xj)= / pjG{x — Xj) dp{x). 


Therefore 


ifi dp = 0 , 


which implies that can be uniformly approximated on S by elements of the 
subspace U, due to a consequence of the Hahn-Banach theorem. Accordingly, 
there is a finite set of points {zj}d^i in 'Mf\B and complex numbers Cj such that 
the function 

f>2{x) :=^2cjG{x-Zj) 

1=1 
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approximates the function Lpi uniformly in S: 

\\‘f2 - V5l||cO(S) < 

The lemma then follows by setting 

.]' j' 

(P2{x)'-^CjG^x-Zj) + ^CjG{x +Zj) =: ^ CjG{x-Zj) 
i=i i=i j=-J' 

where cq = 0, C-j = cj and Z-j = zj. Indeed, since S' is a symmetric set, we have 
that, for all x £ S, 

, , , , , , (pi{x)+lfi{-x) (p2{x)-lfi{x) , ip2{-x) - (Pl(-X) 

ip2{x) - ipi{x) = ip2{x} - ^ + -- -2 : 

which implies the desired estimate 

\W2 - ‘/5i||co(S) ^ ~ ‘^lllc'O(S) + 2^\V2 - </5l||co(S) < ^ ■ 

Notice that we have used the identity </5i(x) = Lpi{—x) □ 


To complete the proof of the lemma, notice that the even complex-valued func¬ 
tion ip 2 constructed in Lemma |4.1| satisfies 

(4.8) A(^2 + ip2=^ 


in the ball B, whose interior contains S. Let us take spherical coordinates {r,9,ip) 
in the ball B. Expanding the function ip 2 (with respect to the angular variables) 
in a series of spherical harmonics and using Eq. (4.8), we immediately obtain that 
Lp 2 can be written in the ball as a Fourier-Bessel series of the form 


OO I 

^2 = EE cimji{r)Yi^{0,(p) . 

l—O m— — l 


Since ip 2 is even, we have that cim = 0 for all odd 1. As before, ji denotes a spherical 
Bessel function. 

Since the above series converges in L^{B), for any (5 > 0 there is an integer Iq 
such that the finite sum 

E (^irnji{r)Yi^{e,(p) 

l—O rn——l 


approximates the function (p 2 in an sense: 

(4.9) llyj-(^2 ||l2(b) < (5. 

By the properties of spherical Bessel functions, the complex-valued function ip is 
smooth in and satisfies the equation 


(4.10) Ap + p = {) 

in the whole space. 

Given any smaller ball B', properly contained in B and in turn containing the 
set S, standard elliptic estimates allow us to pass from the bound (4.9) to a 
uniform estimate 


\W ~ '^2\\c°(B') < C5 . 
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From this inequality and the bounds ( |4.5[ ) and (4.7) we infer 
(4.11) \\ip-ip'\\co(s) <C5. 

Moreover, since Lp' also satisfies the Helmholtz equation in a neighborhood of the 


compact set S (cf. Eq. (4.21), standard elliptic estimates again imply that the 
unifori 

(4.12) 


uniform estimate (4.11) can be promoted to the bound 

llv^- <^'llci(S) < C5. 


Finally, since the link L is a union of components of the the nodal set of ip’ 
and satisfies the transversality condition (|4.1|), the estimate (4.12) and a direct 


application of Thom’s isotopy theorem [U Theorem 20.2] imply that there is a dif- 
feomorphism of such that <l>i(L) is a union of components of the zero set 
(/?“^(0). Moreover, the diffeomorphism is C^-close to the identity and different 
from the identity just in a small neighborhood of L, so we can safely assume that 
$i(L) is contained in B. The structural stability of the link <i>i(L) for the func¬ 
tion p also follows from Thom’s isotopy theorem and the fact that p satisfies the 
transversality condition 

rank(VRev5(T), VIm(/3(a;)) = 2 


for all X € <i)i(L). This last equation is a consequence of the C^-estimate (4.12), 


the fact that the function p' satisfies the transversality estimate (4.11 by definition. 


and the fact that transversality is an open property under C^-small perturbations. 
The lemma then follows. 


5. A REMARK ABOUT THE HIGHER DIMENSIONAL COUNTERPART 


Following Berry, we have considered the construction of a complex-valued eigen¬ 
function (or two real-valued eigenfunctions) of the harmonic oscillator in three di¬ 
mensions with a prescribed nodal set of codimension 2 (that is, a link). It is worth 
mentioning that essentially the same argument enables us to construct n eigenfunc¬ 
tions of the harmonic oscillator in M"’* with a prescribed nodal set of codimension 
n. 


However, a technical condition makes the statement considerably more involved 
in the general case. This condition is associated with the requirement that the 
level set be structurally stable. In the situation covered by the main theorem, the 
structural stability follows from the important equation (4.1), which plays a crucial 
role in the proof. The higher dimensional analog of that relation would then be the 
requirement that 


(5.1) rank(Vi/>i(a;),..., Vi/^„(a;)) = n, 

for all X in the prescribed codimension-n nodal set in where ipi, ... ,ipn would 
be real-valued eigenfunctions of the harmonic oscillator. For this condition to hold, 
a topological obstruction is that the normal bundle of the set L that we want 
to prescribe in the nodal set must be trivial. Geometrically, this is equivalent to 
the assertion that a small tubular neighborhood of the submanifold L must be 
diffeomorphic to L x K". 

Hence we are led to the following result. Since a link always has trivial normal 
bundle m, the main theorem corresponds exactly to the case d = 3 and n = 2. 
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Theorem 5.1. Let L be a finite disjoint union of codimension-n compact subman¬ 
ifolds o /with trivial normal bundle, and d ^ 3. If n = 1, we also assume that 
L is connected. Then for any large enough eigenvalue A of the harmonic oscillator 
in there are n real-valued eigenfunctions ipi,... ,tjjn with eigenvalue A and a 
diffeomorphism $ of such that $(T) is the union of connected eomponents of 
the joint nodal set H • • • n^“^(0). Furthermore, $(T) is structurally stable. 


We note that when d ^ 3 the eigenfunctions of the harmonic oscillator are given 
by 

-ipkim = e~'^r‘L'^i^ ^ {r'^)Yim{uj), 

and the corresponding eigenvalues are X^i = 4fc + 2^ + d. Here ui := x/r is a point 
in and Yim are the spherical harmonics in with frequency l{l + d — 2), 
with I ^ 0, m = (mi,..., md- 2 ) and 


|mi| ^ m 2 ^ < md-2 < I ■ 

Using again Hilb’s asymptotic formula for the Laguerre polynomials 


Theorem 


8.22.4], we get the asymptotic expressions of Lemma 2.1 generalized to any d ^ 3: 

1pklm{x) = Ail [jf (xAfelO + Yi^{uj ), 

Vf)klm{x) = y/^lAii [(jf)'(\/A^r) + Yim{0,(l)) Cr 

+ <, [If (vT^r) + 0{k -»”"<¥.-))] , 

where jf denote the hyperspherical Bessel functions, which satisfy the radial Helmholtz 
equation in and the constants Aii are given by 


Ad — 
■^ki ■— 


1 f y/Xki A ^r(fc + z+|) 


r(f) 


fc! 


The proof of Lemma |2.2| goes exactly as in Section when the codimension is 
n ^ 2. In this case, we use the Cauchy-Kolaweski theorem with data on hypersur¬ 
faces (1 ^ ^ n), intersecting at a given component La transversally, to define 

functions u{,. When n = 1, however, the construction of the real-valued function u 
(analogous to the functions considered in Section where we are dropping the 
sub- and superscripts because L is now connected and n = 1) cannot be performed 
using the Cauchy-Kowalweski theorem because otherwise u would be defined in a 
small neighborhood U of L. As the complement of U would not be connected, the 
proof of Lemma |4.1| does not carry over to this case. To circumvent this difficulty, 
we take Lf to be the precompact domain bounded by L and define u in s 17 as 
the first Dirichlet eigenfunction of the Laplacian in this domain, where s is a scale 
factor chosen so that the first Dirichlet eigenvalue oi sU is 1. The details are as 
in 0 Appendix A]. The rest of the proof remains essentially unchanged. 

In particular, in three dimensions the general result yields not only Theorem |1. 1 1 
but also the following 


Corollary 5.2. Let L be a compact surface in Then for any large enough 
eigenvalue A of the harmonic oscillator in there is a real-valued eigenfunction 
if of the harmonic oscillator with energy A and a diffeomorphism $ of such 
that ^(L) is a connected component of the nodal set if~^{0). Furthermore, $(L) is 
structurally stable. 
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